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ABSTRACT  OF  THE  DISSERTATION 
Discovery  of  a  Nonpropagating  Hydrodynamic  Soliton 

by 

Jon-Rn  Vn 

Doctor  of  Philosophy  in  Physics 
University  of  California,  Los  Angels,  1985 
Professor  Isadore  Rudnick,  Chairman 

A  new  type  water-surf ace-wave  soliton,  nonpropagating 
hydrodynamic  soliton,  has  been  discovered^ .  This  is  a 
self-trapped,  highly  localized  and  stationary  transverse 
water-surface-wave  excitation  which  appears  in  a  rectangular 
cross  section  resonator  continuously  excited  parametrically  by 
vertical  oscillation.  ) (Fig.  1). 

An  experiment  designed  to  create  the  soliton  and  the 
features  of  the  soliton  are  described.  The  profiles  of  the 
soliton  for  various  frequencies  and  amplitudes  are  measured.  It 
is  found  that  the  profile  of  the  soliton  is  accurately  given  by  a 
hyperbolic  secant  function.  The  stability  region  that  is  the 
amplitude  and  frequency  range  of  the  drive  in  which  individual 
soliton  can  be  created  without  hysteresis  is  experimentally 
determined. 

The  interaction  between  two  solitons  has  been  investigated. 
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It  it  observed  that  an  attractive  force  exists  between  two 
solitons  with  the  sasie  polarity  and  a  repulsive  force  exists 
between  two  solitons  with  the  opposite  polarity.  Under  a  certain 
range  of  drive  frequencies  and  amplitudea,  two  solitons  with  the 
sane  polarity  slowly  oscillate  through  each  other.  The 
oscillation  frequency  has  been  neasured. 

The  self-trapping  aechanisn  of  the  soliton  is  studied.  It 
is  found  that  the  necessary  condition  for  creating  the  soliton  is 
that  our  nonlinear  dispersive  system  has  an  amplitude  dependent 
resonant  frequency  of  the  transverse  mode  such  that  the  higher 
the  amplitude  the  lower  the  resonant  frequency.  The  finite 
amplitude  tuning  curves  of  the  rectangular  cross  section 
water-surface-wave  have  been  determined  for  varioua  liquid 
deptha.  The  experimental  results  show  that  the  direction  of  the 
bending  of  the  tuning  curve  depends  on  the  liquid  depth.  For  a 
deep  liquid  it  bends  toward  lower  frequencies,  and  for  a  shallow 
liquid  it  bends  toward  higher  frequencies.  The  so-called 
nonlinear  coefficient  is  determined  from  these  tuning  curves. 
The  cross-over  from  right  leaning  to  left  leaning  was  found  to  be 

in  good  agreement  with  the  theoretical  value  calculated  from  the 

(2) 

formula  developed  independently  by  J. Niles  ,  and  by  A.Larraxa 
and  S.  Putterman^. 

The  soliton  free  decay  rate  is  also  measured.  It  is  found 
that  the  soliton  initially  decays  much  faster  than  the 
corresponding  linear  transverse  mode,  but  has  the  same  rate  in 


CHAPTER  I 


INTRODUCTION 


''....it  hat  usually  been  thought 
that  in  a  relativistic  field  theory, 
in  order  to  have  stationary  bound 
states,  quantum  mechanics  nust  be 
crucial.  As  we  shall  see.  this  turns 
out  not  to  be  the  ease.  In  a 
nonlinear  field  theory,  with  an 
appropriate  amount  of  nonlinearity, 
stable  bound  states  can  exist  on  a 
classical,  as  well  as  quantum 
mechanical,  level.  Such  bound  states 
are  called  solitons.'' 


AN  HISTORICAL  REVIEW  OF  SOLITONS  IN  A  FLUID 

The  concept  of  a  solitary  wave  was  first  introduced  by  a 


Victorian  engineer  John  Scott  Rnatell  in  1844  in  his  historical 
paper"  Report  on  Waves"  (appendix  A).  He  nsed  the  word 
"solitary"  to  characterize  the  robust,  localized,  and  well 
defined  "heap"  of  water  which  propagates  along  a  narrow  canal 
without  change  of  fora  and  diminution  of  speed^.  Russell's 
observation  appeared  to  contradict  the  famous  contemporary 
scientist  Airy's  shallow  water  theory  which  predicts  that  a  wave 
of  finite  amplitude  cannot  propagate  without  change  of  form^ . 

The  Airy  paradox  wasn't  totally  resolved  until  1895  when 

Korteweg  and  deVries  developed  an  equation  for  a  water  wave  which 

(7) 

includes  both  nonlinear  and  dispersive  effects  .  This  equation 
has  permanent  solitary  wave  solutions (appendix  B) .  It  is  now 
well  known  as  the  KdV  equation,  and  provides  a  simple  analytical 
foundation  for  the  study  of  solitary  waves. 

A  new  is^etus  to  this  problem  was  generated  by  Fermi,  Pasta, 

(g) 

and  Ulan's  work  in  the  early  fifties  .  They  used  the  first 
large  electronic  computer.  Maniac  I,  to  investigate  the  time 
evolution  of  the  energy  of  a  weakly  and  nonlinearly  coupled 
system  of  sixty  four  harmonic  oscillators.  They  started  with  the 
case  when  the  energy  was  initially  concentrating  in  one  of  the 
modes.  To  their  great  surprise,  the  classical  equipartition 
theorem,  describing  the  thermal  equilibrium  ,  turned  out  to  be 
incorrect.  After  some  tens  of  thousands  of  cycles,  the  energy 
invariably  returned  almost  completely  to  the  original  mode.  This 


is  veil  known  now  «s  F.P.U.  Recurrence  ,  and  is  evidence  that 
the  development  of  snoh  a  collective  node  is  a  general 
phenomenon.  It  occurs  in  the  other  fields  of  the  physics, 
outside  of  hydrodynamics (appendix  C). 

In  1 965  N.  J.  Zabusky  and  M.  D.  Eruskal  found  from  their 
numerical  simulation  of  the  nonlinear  discrete  mass  string 
problem  that  an  intially  sinusoidal  profile  would  decompose  into 
a  series  of  interacting  solitary  wavea  and  each  one  of  them 
retained  its  identity  after  each  interaction1  .  It  is  the  first 
time  that  the  term  soliton  was  used  to  characterize  this  kind 
of  solitary  wave  which  has  the  property  of  emerging  from  the 
collision  having  the  same  shape  and  velocity  with  which  it 
entered  -  a  remarkable  resemblance  to  the  elastic  interaction  of 
particles. 

Although  the  term  soliton  was  originally  applied  only  to  the 
solitary  waves  of  the  KdV  equation,  there  are  now  several 
nonlinear  partial  differential  equations  known  to  exhibit  similar 
effects.  Among  them,  especially  for  fluid  solitons,  the 
Nonlinear  Schrodinger  Equation  (NLS)  is  a  popular  one^*^.  in 
1968  Zakharov  predicted  that  the  time  evolution  of  weakly 
nonlinear  deep-water-waves  obeys  NLS  and  three  years  later 
Zakharov  and  Shaba t  found  the  exact  solution  for  the 
equation^*).  In  197S  Tuen  and  Lake^^  verified  the  prediction, 
and  observed  the  so-oalled  envelope  soliton  in  a  narrow  water 


tank (appendix  D).  The  envelope  aoliton  can  be  pictured  aa  a 
quickly  oscillating  wave  cut  off  by  a  saoothly  modulating 
envelope  (Fig. 2).  Until  recently  the  KdV  and  envelope  solitons 
were  the  only  known  solitary  waves  in  fluids.  Both  of  these 
solitons  have  to  propagate  at  a  non-zero  group  velocity 
determined  by  their  amplitudes  and  shapes.  They  are  also  limited 
to  motion  in  one  dimension. 

Recent  theoretical  efforts,  as  far  aa  fluid  mechanics  is 

concerned,  concentrate  on  extending  the  one  dimensional  classical 

(13) 

soliton  solution  to  higher  dimensions  . 

B.  AN  INTRODUCTION  TO  TUB  NEWLY  DISCOVERER  NOWRQPAGATING 
HYDRODYNAMIC  SOLITON 

The  nonpropagating  soliton  was  first  observed  in  Professor 
I.  Rudnick's  laboratory  in  the  summer  of  1983  in  a  fluid  trough 
driven  sinusoidally  and  vertically  from  below  so  that  all  the 
fluid  is  caused  to  accelerate  up  and  down  equally.  A  localized 
surface  pulse  could  be  observed.  The  motion  of  the  pulse 
appeared  to  be  so  stable  and  robust  that  it  was  immediately 
thought  to  be  a  soliton.  Furthermore,  since  it  could  appear  and 
be  stationary  anywhere  along  the  length  of  the  trough  and  sloshes 
back  and  forth  across  the  width  of  the  trough,  more  than  one 
dimension  is  involved.  It  is  fundamentally  different  from  the 


other  two  types  of  soliton  and  is  called  a  nonpropagating 
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To  understand  the  physics  of  a  nonpropagating  soliton,  it  is 


itial  to  review  waveguide  theory.  Suppose  the  trough  is  very 


long  so  that  it  can  be  regarded  aa  a  waveguide  with  length 
coordinate  z  and  width  coordinate  y  (Fig. 3).  The  depth  of 
liquid,  d.  controls  the  surface  wave  phase  velocity  ,c,  given  by 


c  *  (g/k  +pk/p)tanh(kd) 


(  1  ) 


where  k  is  the  wave  number,  g  is  the  acceleration  of  gravity,  p 
is  the  surface  tention  and  p  is  the  density  of  fluid.  The 
surface  displacement,  z,  is  given  by 


■  a  cos(k^y)exp[i(kxx  *  2nvt)] 


(  2  ) 


where  v  is  frequency,  k^W  ”  qn,  W  is  the  width  of  the  trough,  q 
0,1,2,3...  and 


kx  "  <  k*  -  V  >’  -  (2ji/c)  (  v*  -  vo;toff  )’ 


cutoff  2W 


(  3  ) 


(  4  ) 


The  case  q  *  0  is  the  plane  wave  mode,  q  ■  1  has  one  velocity 


9 


component,  v  ,  antinode  between  tbe  side  wells,  q  *  2  has  two, 
etc.  The  energy  of  the  wave  propagates  with  the  gronp  velocity, 
▼  ,  which  can  be  calculated  by 


c 

g 


v 

V 


cutoff 


(  5  ) 


If  v  >  vcnt0ff»  and  Cg  are  real,  energy  propagates  down  the 

.1  %  2 

waveguide  in  a  angle  of  0  *  cos  [1  -(v^  ogj/v)  1  with  the  x 
direction(Fig.  4).  If  v  <  vent0f£  then  kx  and  Cg  are  imaginary. 
The  amplitude  of  this  particular  mode  decays  exponentially  down 
the  waveguide,  and  energy  will  not  propagate,  but  is  instead 
totally  reflected  back.  This  is  a  so-called  evanescent  wave. 
For  a  finite  waveguide  of  length  L  like  our  resonator. 


*  ■  a  cos(kix)cos(ky.y)cos(2jtvt )  (  6  ) 

where  k^L  “  pn  and  p  *  0,1,2,3....  Resonant  frequencies  are 
given  by 

x 

v  „  ■  <c/2)[(p/W)#  +  (q/L>*]\  (  7  ) 

P*R 

The  (0,1)  mode  is  the  equivalent  of  being  the  first  cutoff  for 


the  infinite  waveguide 


OBSERVATION  OF  A  NONPROPAGATING  SOL I TON 


' '  Ton  sake  experiment*  and  I  make 

theories.  Do  yon  know  the  difference? 
A  theory  is  something  nobody  believes 
except  the  person  who  made  it,  while 
an  experiment  is  something  everybody 
believes  except  the  person  who  made 
it." 

A.  Einstein^**^ 


The  apparatus  for  getting  a  nonpropagating  soliton  is  quite 
simple.  Te  use  a  Plexiglas  (or  glass)  channel  38  cm  long  and  2.6 
cm  wide  filled  with  still  water  (  unless  otherwise  specified,  the 
working  liquid  is  still  water)  to  a  depth  of  2  cm.  Several  drops 
of  the  wetting  agent  Kodak  Photo-Flo  are  added  to  the  water  to 
minimize  surface  pinning  effects  at  the  walls.  Waves  are 
generated  by  placing  the  trough  horizontally  on  a  18" 
loudspeaker  (or  a  vibration  exciter,  4809,  Bruel  Ijaer)  whose 
cone  is  driven  at  a  frequency  2v  in  the  vertical  direction.  As 


•.own  by  Eq. (  7  ) ,  there  could  be  (p,q)  modes  appearing  in  the 
trough. 


It  is  found  thst  with  this  psrsaetric  excitation  all 
low-lying  (p.O)  modes  can  be  observed.  When  we  parametrically 
excite  the  (0.1)  mode,  however,  there  is  some  degree  of 
competition  between  various  modes  at  the  initial  stage,  then  the 
sloshing  motion  in  y  direction  corresponding  to  (0,1)  mode 
developes  and  appears  to  be  unstable.  Finally  the  waves  merge 
and  focus  themselves  to  one  or  more  excitations  highly  localized 
in  x-direction,  sloshing  in  y-direction.  Importance  is  attached 
to  this  self-focussing  behavior.  It  is  surprising  and  intriguing 
that  while  all  parts  of  the  trough  are  oscillated  with  equal 
amplitude,  the  vigorous  wave  motion  occurs  in  a  space  of  a  few 
centimeters  long  while  the  rest  of  the  water  remains  tranquil. 


The  phase  shift  along  the  x-direction  given  by  Eq.(  2  )  for 


vcutoff  “  5,1  Mnd  v  “  5,4  i#  90  in  *  di9ttnce  of  6  c** 


Tet,  within  our  resolution  of  3  ,  we  find  no  measurable  phase 
difference  between  any  two  points  of  the  disturbance.  It  is 
believed  that  the  absence  of  a  phase  difference  is  of  fundamental 
importance.  The  wave  is  not  simply  a  waveguide  mode  occuring 
slighly  above  the  cutoff  frequency  but  a  new  type  of 
nonpropagating  solitary  wave.  As  we  proceed  we  shall  see  that  it 
has  properties  usually  associated  with  solitons. 


Figure  5  and  6  are  computer  generated  profiles  of  the 


soliton  based  on  aeasureaents  taken  froa  photographs,  when  it  is 
respectively  at  its  peaks  on  the  far  side  of  the  trough  and  on 
the  near  side  of  the  trough.  The  wave  frequency  here  is  5.1  hz 
and  the  peak  aaplitude  is  2.1  ca. 

Figs  7  -  10  contain  plots  of  the  heights  z  of  the  solitons 
as  functions  of  z  for  various  frequencies  and  the  same  aaplitudes 
of  the  drive,  when  they  are  at  their  peaks  on  the  vails.  The 
points  are  ezperiaental  data  and  the  curves  are  least  square 
fits.  Froa  the  plots,  it  is  evident  that  our  solitons  have  the 
following  features. 

The  hyperbolic  secant,  a  function  vhich  characterizes  soae 
solitons,  fits  our  soliton  profile  in  z  direction  very  veil.  The 
aaplitude  of  the  nonpropagating  soliton  can  be  as  high,  or  higher 
than  the  water  depth. 

The  characteristic  length,  vhich  is  defined  as  the  characteristic 
decay  length  at  its  vings,  and  the  aaplitude  of  the  soliton  are 
deterained  by  the  frequency  of  the  ezeitation  for  the  saae 
aaplitude  driven  systea.  They  are  interrelated  in  such  a  way 
that  the  higher  the  aaplitude  the  shorter  the  characteristic 
length. 

The  frequency  of  the  sloshing  notion  of  the  nonpropagating 
soliton  is  lover  than  v  given  by  Eq.  4. 

The  dependence  of  z  on  y  vhen  the  soliton  is  at  its  peak  is 
found  to  be  far  froa  a  cosine  function  which  is  linear  solution 


to  the  Laplace  equation.  A  beat-fit  function  is  found  to  be  an 
exponential  function  (Fig.  11) . 


It  is  found  experimentally  that  one  of  the  important  feature 
of  our  soliton  is  that  the  velocity  along  the  x-direction  of  the 
localized  disturbance  is  a  free  parameter.  It  can  be  zero  if  the 
trough  is  accurately  horizontal,  which  means  the  soliton 
maintains  its  initial  x-position  in  the  trough  for  the  order  of 
an  hour  or  more.  Ve  hypothesize  that  this  stability  is  due  to 
pinning  effects  at  the  wall.  The  soliton.  however,  can  be  moved 
in  various  ways.  One  such  way  is  to  tilt  the  trough.  For  a 
slope  of  0.05,  for  instance,  the  soliton  slowly  moves  to  the 
shallow  end  with  an  average  speed  of  0.05  cm/sec.  It  moves  in 
response  to  gentle  jets  of  air.  It  cun  also  be  nudged  by  rods  of 
soft  sponge  plastic.  Unwanted  solitons  are  killed  by  stabbing 
them  with  such  rods. 

The  stability  range,  defined  as  the  range  of  the  amplitudes 
and  frequencies  of  the  drive  in  which  individual  solitons  are 
observed  without  hysteresis,  were  measured  and  found  to  be  as 
shown  by  the  dashed  lines  of  Fig.  12.  The  eight  full-line  curves 
are  equal-wave-response  curves.  The  number  on  each  is  the  peak 
height  in  centimeters  of  the  soliton  above  the  equilibrium  level 
of  the  water.  Note  that  at  0.7  and  2.1  cm  the  curves  are  so 
short  that  they  are  only  represented  by  points.  From  the  plot  we 
see  that  for  the  same  amplitude  of  drive,  the  lower  the  frequency 
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experiment*!  date,  and  the  aolid  «me  i*  x  “  2.1*ech(0.89x)  cm 


the  higher  the  amplitude  of  the  soliton.  Outside  the  dashed 
boundary,  changes  in  drive  amplitude  and/or  frequency  are 
generally  hysteretic.  If  it  exists,  the  soliton  may  be  unstable, 
and  there  may  be  a  mixture  of  modes.  At  a  drive  amplitude 
exceeding  0.081  cm,  charateristically  a  competing  mode  rather 
than  a  soliton  appears.  As  mentioned  before,  there  is  generally 
some  degree  of  competition  between  the  various  modes  of  the 
trough,  especially  after  the  drive  is  just  turned  on.  However, 
if  the  amplitude  and  frequency  of  the  drive  are  appropriately 
set,  energy  will  finally  be  sucked  into  a  soliton  and  it  becomes 
stable.  Several  ways  are  devised  to  selectively  encourage  the 
appearance  of  the  soliton  to  the  exclusion  of  competitive  modes. 
One  obvious  way  is  to  choose  the  trongh  dimensions  which 
separate,  in  frequency,  the  unwanted  modes  from  the  soliton  mode. 
Another  way  is  to  put  sponge  dampers  at  both  ends  of  the  trough 
so  that  the  longitudinal  modes  are  discouraged.  Another  very 
successful  procedure  is  to  nucleate  the  soliton  mode  by  producing 
a  local  disturbance  which  is  compatible  with  it.  Thus  sloshing 
motion  across  the  width  can  be  produced  by  rocking  the  resonator 
or  with  the  help  of  a  hand-held  paddle.  This  can  be  very 
effective  if  it  is  done  when  the  drive  is  just  turned  on.  The 
self-focusing  process  described  before  is  very  striking  if  the 
soliton  is  generated  by  rocking  the  trough. 

Since  the  frequency  of  the  wave  is  half  that  of  the  drive, 
the  phase  of  the  drive  is  repeated  every  180*of  the  wave. 


Consequently  solitons  which  ere  0  or  180  apart  in  phase  are 

coapatible  and  equally  driven.  Thus  the  only  steady  state 

• 

soliton  pairs  are  in-phase  solitons  of  the  saae  polarity  or  180 
out-of-phase  solitons  of  the  opposite  polarity.  An  opposite 
polarity  soliton-pair  can  be  generated  by  sharply  twisting  the 
trough  about  a  central  vertical  axis,  so  that  the  sloshing  in  the 
left  part  has  the  opposite  direction  of  that  in  the  right  part. 
Alternatively  an  inverted  T-shaped  paddle  iaaersed  in  the  water 
can  be  twisted  to  achieve  the  saae  result. 

The  interaction  of  the  solitons  was  studied.  Two  solitons 

of  the  saae  polarity  attract  each  other,  but  only  weakly  if  the 

distance  between  then  is  significantly  larger  (say  a  factor  of  3) 

than  their  characteristic  length.  Two  solitons  which  start  out 

20  ca  apart  center  to  center,  for  exaaple,  take  about  15  ain  to 

reach  a  separation  where  they  strongly  overlap.  When  this  happen 

their  attractive  speed  greatly  increases.  If  the  frequency  is 

substantially  lower  than  the  saall  aaplitude  v#  x,  the  two 

» 

solitons  coabine;  the  end  result  is  a  single  soliton  having  the 

aaplitude  of  each  of  the  initial  ones.  If  the  frequency  is 

closer  to  v0  x,  the  solitons  oscillate  through  each  other.  Fig. 
» 

13  is  a  scheaatic  of  the  stages  of  the  oscillation.  In  the  first 
stage.  Fig. 13(a),  there  is  significant  overlap.  In  Fig. 13(b)  the 
overlap  is  sufficient  to  create  an  aaplitude  at  the  center  which 
is  coaparable  with  the  aaplitudes  of  solitons  1  and  2.  In 
Fig. 13(c)  the  solitons  coapletely  overlap.  In  Fig. 13(d) 


apparently  sol itons  1  and  2  have  passed  through  each  other  and 
exchanged  places  while  in  Fig. 13(e)  the  configuration  is  the 
original  one.  Left  undisturbed  this  sequence  repeats 
indefinitely.  The  longest  period  of  tine  in  which  this 
phencaenon  was  observed  is  36  hours.  An  interruption  of  the 
oscillation  has  never  been  seen  except  by  accident  or  design. 

A  simple  transducer  which  responds  to  wave  heights  is  a  pair 
of  vertical  wire  electrodes  that  dip  into  the  water.  At  a 
constant  ac  voltage  drive,  the  current  due  to  the  conductance  of 
the  water  is  proportional  to  the  wave  height  at  the  electrodes. 
Fig.  14  shows  the  response  of  such  a  transducer  placed  at  the 
midpoint  of  the  oscillating  soliton-pair  of  Fig.  13.  In  the 
lowest  trace,  which  covers  the  complete  sequence  of  Fig. 13,  the 
paper  speed  was  great  enough  so  that  the  individual  oscillations 
are  seen  (v  =  3.32  hz).  The  frequency  of  the  envelope,  which  is 
much  lower  than  the  drive,  is  2v/(156±0.5)  =  0.068  hz.  The  upper 
traces  are  recorded  at  succesively  slower  paper  speeds.  Fig.  15 
is  the  identical  plot  for  the  case  that  the  fluid  is  pure 
alcohol.  The  principal  difference  is  that  the  frequency  of  the 
oscillation  is  a  little  higher  (2v/ (103±0.5 )  =  0.103  hz  )  than 
the  water  case.  This  suggests  the  surface  tension  might  be 
playing  an  important  role  in  the  oscillation.  The  surface 
tension  of  alcohol  is  only  one  third  that  of  water. 

A  pair  of  solitons  of  opposite  polarity  in  close  proximity 


to  each  other  repel  each  other  and  slowly  move  until  they  are 
approximately  12  cm  apart,  and  then  maintain  this  separation 
indefinitly.  Again  the  interruption  of  this  state  has  never  been 
seen  except  by  accident  or  design. 

By  introducing  tracers  (eccospheres)  into  the  water  we  have 
been  able  to  establish  that  the  motion  of  solitons  is  one  in 
which  there  is  a  transport  of  excitation  with  negligible 
transport  of  mass. 

Essentially  identical  phenomena  have  been  observed  in  an 
annular  resonator  72  cm  in  mean  circumference  and  2.2  cm  wide. 
This  is  interpreted  to  mean  that  the  boundary  conditions  in  the 
x-direction  are  irrelevant  to  our  problem.  As  many  as  six 
solitons  have  been  observed  in  the  resonator.  The  polarity  of 
solitons  always  alternates  in  the  steady  state  limit. 


idpoit  of  the  oscillating  solitons  of  Fig. 


CHAPTER  III 


A  NONLINEAR  PHYSICAL  MODEL  OF  THE  NONPROPAGATING  SOLITON 


"  The  sol it on  digs  its  own  hole 
sad  builds  its  own  wells.*' 

I.  Rudnick 


The  experimental  observations  introduced  in  Chapter  II  can 
be  summarized  as  follows. 

(1)  The  system  is  very  nonlinear  as  the  amplitude  of  the  solitoa 
can  be  as  high  as  water  depth. 

(2)  The  profile  of  the  soliton  in  the  x  direction  can  be 
described  by  the  hyperbolic  secant  function,  which  asymptotically 
approaches  exponential  decay  function  at  its  tails. 

(3)  The  frequencies  for  getting  the  nonpropagating  soliton  is  in 
a  frequency  range  lower  than  the  cutoff  frequency  given  by 
Eq.(  4  ).  Moreover  the  higher  the  drive  amplitude,  the  lower  the 
optimal  frequency  for  peak  response. 

(4)  The  nonpropagating  soliton  has  the  tendency  to  move  to  the 
shallow  end  where  the  cutoff  frequency  is  lower. 

What  kind  of  the  physical  model  can  we  use  to  describe  the 
nonropagating  solitoa  based  on  these  experimental  facts!  A 
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reasonable  model  that  would  eome  up  is  that  the  cutoff  frequency 
of  the  system  becomes  amplitude  dependent  as  the  system  appears 
to  be  nonlinear.  The  fact  that  the  drive  frequency  for  getting 
the  nonpropagating  soliton  is  below  the  cutoff  frequency  suggests 
the  model  should  be  such  that  the  higher  the  amplitude  of  the 
waves  the  lower  the  cutoff  frequency.  The  simplest  system  which 
has  this  kind  of  property  is  an  undamped  free  oscillating 
pendulum.  When  friction  is  neglected,  the  differential  equation 
governing  the  motion  of  the  pendulum  is  given  by 


*  +  uc*sinz  *  0,  (  8  ) 

where  *  g/L;  the  angle  z  designates  the  deviation  from  the 
vertical  equilibrium  position  (Fig.  16).  This  is  a  nonlinear 
differential  equation  as  we  know  the  Taylor  Series  of  sinx  is 
given  by 


sinz  ■  z  -  z*/3!  +  z*/5!  -  z*/7!  + 


(  9  ) 


Only  in  the  first  order  approzination,  does  Eq.  (  8  )  reduce  to 
the  linear  differential  equation: 


Then  the  resonant  frequency  of  the  system  is  independent  of  the 


angle, 
justified 
Eq.  (  9 


This  is 
order)  of 


where  <■» 
equation: 


For  large  angles  the  linear  approximation  is  not 
If  the  second  term  of  the  expansion  is  also  included, 
becomes 


x  +  «#*x  -  («#*/6)x*  *  0.  (11) 

the  Duffing  equation*15* .  The  solution  (up  to  third 
the  nonlinear  differential  equation  can  be  written  as 

x  ■  asin(ut  +y)  +  t*a*/(192u*)a>sin(3ut  +  3y),  (12) 

is  amplitude  dependent  and  is  given  by  the  following 


«*(a) 


1  /  » 
“•  “  /*)»  • 


(13) 


Now  let  us  use  this  model  to  understand  the  self-trapping 
effects  of  the  nonpropagating  soliton.  At  the  core  of  the 


sol itoa  the  peek  anplitude  st  the  wall  is  so  high  that  the 
corresponding  entoff  frequency  there  is  significantly  lower  than 
the  snail  amplitude  entoff  frequency  given  by  Eq.  (  4  ) 
according  to  our  nodel.  However  the  linear  theory  is  still  valid 
at  both  wings  of  the  nonpropagating  soliton  since  the  aaplitudes 
there  are  pretty  snail.  If  the  systen  is  driven  at  the  frequency 
higher  than  the  anplitude  dependent  cutoff  frequency  at  the  core 
and  lower  than  the  snail  anplitude  frequency  at  the  wings,  the 
energy,  which  propagates  out  of  the  core,  gets  totally  reflected 
at  the  wings,  and  the  wave  becones  evanescent.  This 
self-trapping  is  reniniscent  of  polaron  behavior  in  crystals,  and 
leads  to  the  soliton  being  called  a  hydrodynamic  polaron. 

Consider  Fig.  7.  The  soliton  profile  in  the  x  direction  is 
described  by  z  *  2.1sech(x/1.12) cn.  At  the  core  the  peak 

anplitude  at  the  wall  is  2.1  on  ,  consequently  the  cutoff 
frequency  is  lower  than  that  at  the  wings.  At  the  tails  of  the 
hyperbolic  secant  function  the  decay  length,  1.12  cn,  is  equal  to 
the  low  anplitude  waveguide  evanescent  length,  and  this  is  just 
the  nagnitude  of  1/k^.  Thus  we  can  deternine  the  phase  velocity, 
c,  by  the  following  equation: 

i 

c»»/(k*+k*) 

x  y  i 

-  2nx5.1/ [ (n/2.6) 1  -  (1/1. 12)*]* 

»  37.7  cn/sec. 


Eq.  (  1  )  gives  the  value  33.7  ca/sec.  This  siaple  calculation 
shows  that  our  aodel  is  reasonably  good. 

By  using  this  aodel  it  is  not  difficult  to  understand  the 
self-focusing  effects  observed  at  the  initial  stage  of  the 
process  of  the  soliton  formation.  As  the  drive  aaplitude  passes 
the  threshould  value  of  the  paraaetrie  instability,  the  aaplitude 
of  the  excitation  starts  to  develops.  Wherever  the  aaplitude  is 
higher,  the  aaplitude  dependent  cutoff  frequency  becoaes  lower. 
As  has  already  been  pointed  out  that  the  energy  of  the  excitation 
likes  to  go  to  the  place  where  the  cutoff  frequency  is  lower.  In 
a  snowball-like  process,  the  energy  builds  up  in  a  very  saall 
region,  and  the  soliton  is  foraed.  Once  the  soliton  appears  it 
is  stable  due  to  the  above-mentioned  self-trapping  mechanism. 
Another  fact  which  is  ia^ortant  to  the  stability  of  the  soliton 
is  the  dispersion  of  the  systea.  If  there  is  no  dispersion, 
shock  waves  fora  and  there  is  abundant  energy  dissipation. 


MEASUREMENTS  OF  TUNING  CURVES 


As  pointed  out  earlier  the  necessary  conditions  for  creating 
a  nonpropagating  soliton  is  that  a  nonlinear  dispersive  system 
has  an  amplitude  dependent  cutoff  frequency  such  that  the  higher 
the  amplitude  the  lover  the  cutoff  frequency.  Therefore  in  order 
to  test  our  physical  model  it  is  necessary  to  measure  the 
amplitude  dependent  resonant  frequency  of  (0,1)  mode  of  our 
trough. 

A.  INTRODUCTION 

The  amplitude  dependent  resonant  frequency  of  (0,1)  mode 

(2) 

calculated  by  J.  Miles  ,  and  by  A.  Larraza  and  S. 

(3) 

Putterman  is  given  by  the  following  equation  (  Eq.  (E-24), 
Appendix  E  ) : 

v,(0  -  v,[l-AgY/(128wY)]  (14) 

where  v,  is  determined  by  Eq.  (  4  ),  and  A  is  given  by  Eq.(E-20) 


(Appendix-E)  as 


A  =»  |  (  6T4  “  5T*  +  16  -  T~*  )  (15) 

where  T  =  tanh(^d).  Then  correction  is  Bade  for  surface  tension 
Eq.  (E-20)  becones  Eq.  (E-22)  as 

A  *  y[l+(l-T*) *+2(l+44*)-1(l+T*) *  -J(T*-44*)-1(3-T,),j  (16) 

where  p  is  surface  tension,  and  p*  *  p(^)*/(pg).  Proa  Eq.  (14) 
it  is  obvious  that  the  aaplitude  dependent  frequency  v  can  be 
greater  or  less  than  v0  depending  on  the  sign  of  A,  which  is  in 
turn  uniquely  deternined  by  the  depth  of  the  liquid.  By  putting 
A  -  0  in  the  equation  the  theoretical  crossover  depth  can  be 

calculated.  Experiaentally  v,(£)  can  be  deternined  by  neasuring 
the  tuning  curve  of  the  (0,1)  node  for  various  amplitudes,  (,  of 
the  surface-wave  oscillation  for  different  depths  of  liquid.  By 
using  the  measured  value  of  v#  -  the  peak  response  frequency  at 
small  aaplitude,  v#((;)  -  the  peak  response  frequency  at  high 

aaplitude  ,  £,  when  we  sweep  the  frequency  for  different  drive 

amplitudes,  the  magnitude  as  well  as  the  sign  of  the  nonlinear 
coefficient.  A,  can  be  determined  from  the  Eq.  (14). 


The  resonator  used  to  measure  the  toning  curve  is  a  sealed 
glass  basin  with  the  inside  dimensions  of  L  2.0  x  W  2.6  x  H  7.0 
in  centimeters. 

The  reason  for  not  osing  the  long  channel  (38  cm)  to 
determine  the  toning  curve  is  that  when  the  amplitude  of  the 
surface-wave  gets  quite  high,  solitons  appear  in  the  long 
channel.  The  2cm  length  is  chosen  to  be  smaller  than  the 
characteristic  length  of  the  soliton  so  that  the  soliton  does  not 
appear  and  so  that  the  fundamental  resonant  frequency  of  the 
longitudinal  (1,0)  mode  is  well  above  that  of  the  (0,1)  mode. 
Then  the  tuning  curve  is  that  of  the  transverse  (0.1)  mode,  and 
as  pointed  out  earlier,  this  is  the  equivalent  of  being  at  the 
first  cutoff  frequency  for  a  waveguide  of  infinite  length. 

The  resonator  is  cleaned  with  a  potassium  hydroxide,  ethyl 
alcohol,  solution  followed  by  washing  with  hot  tap  water  and 
ethyl  alcohol.  Then  it  is  filled  with  100%  ethyl  alcohol  up  to 
the  depth  of  d  as  the  working  fluid.  The  resonator  is  put  on  an 
improvised  table  on  rollers  which  is  driven  horizontally  by  a  5 
inches  loudspeaker  in  a  direction  normal  to  the  sides  of  the 
basin  so  that  sloshing  occurs  across  the  width  of  the  basin  (Fig. 
13).  An  accelerometer  (4367,  Bruel  Ejaer)  is  used  to  measure 
the  displacement  of  the  oscillating  table.  It  is  found  that  the 
displacement  is  quite  sinusoidal  —  the  second  hsrmonic  component 
is  at  least  40  db  lower  than  the  fundamental. 
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C.  ELECTRONICS 


The  block  diagram  of  the  eletronics  is  shown  in  Fig.  16. 
The  oscillator  (Wavetek)  is  slowly  swept  over  the  range  of  3  -  7 
hz  by  using  an  external  ramp  voltage  signal.  The  output  of  the 
oscillator  drives  the  loudspeaker  after  going  through  a  power 
amplifier.  In  order  to  keep  the  displacement  amplitude  of  the 
oscillating  table  constant  when  the  frequency  is  changed,  a 
feedback  loop  is  used.  The  basic  principle  is  as  follows.  An  ac 
voltage  signal  from  the  accelerometer  proportional  to  the 
oscillating  table  displacement  is  sent  to  the  charg  amplifier 
(2635,  Bruel  4  Kjaer)  then  is  rectified  and  filtered  by  an  ac  - 
dc  converter.  A  dc  differential  amplifier  is  used,  its  negative 
input  is  connected  to  the  output  of  the  ac  -  dc  converter  and  its 
positive  input  is  controlled  by  a  dc  regulated  power  suply.  Thus 
the  amplified  negative  dc  signal  from  the  output  of  the  dc 
differential  amplifier  is  related  to  the  displacement  amplitude 
of  the  oscillating  table  in  such  a  way  that  the  higher  the 
displacement  amplitude  the  more  negative  the  dc  output.  This 
negative  dc  voltage  is  sent  to  the  VCA  input  of  the  Wavetek  to 
decrease  the  amplitude  of  the  ac  output  of  Wavetek.  The  result 
of  this  simple  feedback  loop  is  quite  good;  it  makes  the  driving 
system  which  originally  has  a  peak  in  response  at  about  5  hz  to 
become  a  constant  amplitude  of  displacement  driving  source. 


□ 


The  transducer  for  measuring  the  displacement  of  the 
surface-wave  is  a  pair  of  vertical  wire  electrodes  that  is 
immersed  into  the  liquid.  A  built-in  low  output  impedance 
(600JI)  ac  voltage  (1  v,  1000  hz)  from  lock-in  amplifier  (HR  -  8, 
Princeton  Applied  Search)  is  applied  on  the  transducer,  which  is 
in  series  with  a  100  J1  resister.  Since  the  resistance  of  the 
resister  is  negligible  comparing  with  the  impedance  of  the  liquid 
between  the  electrodes  (about  100  kA),  the  voltage  drop  across 
the  resister  is  proportional  to  the  conductance  of  the  liquid 
between  the  pair  of  electrodes.  When  the  height  of  liquid 
oscillates  at  about  5  hz  as  the  surface-wave  does,  the  ac 
conductance  of  the  liquid  at  the  transducer  is  modulated.  The 
lock-in  amplifier  picks  up  the  slow  modulation  signal,  amplifies 
it  and  sends  it  to  the  ac  -  dc  converter.  An  x-y  recorder  is 
used  to  plot  out  the  tuning  curves. 


Typical  frequency  sweeps  for  several  different  drive 
amplitudes  for  each  of  six  different  depths,  d,  are  shown  in 
Pigs.  19  -  25.  It  is  evident  as  shown  in  Fig.  19  that  when  d  = 
0.5  cm  the  peak  frequency  increases  as  the  amplitude  of 
water-surface-wave  becomes  bigger.  vs  is  4.22  hz  from  Eq. (  4  ) 
in  good  agreement  with  the  observed  4.23  hz.  As  the  amplitude, 
goes  from  0.1  cm  to  0.4  cm  in  steps  of  0.1  cm,  the  peak 
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iplitnde  dependent  toning  curves  of  the  ease  resonator 


iplitnde  dependent 


frequency  shifts  up  from  4.23  hz  to  4.42  hz.  This  is  s  typical 
right  bending  tuning  curve. 

Figs.  20  -  24  are  respectively  for  the  case  of  d  *  1.0. 
1.3.  1.7,  2.0,  and  3.0  cm.  The  coawon  feature  is  that  the  peak 
frequencies  all  shift  down  as  the  amplitudes  of  the  responses  go 
higher  and  higher;  all  belong  to  left  bending  tuning  curves. 
Fig.  20  (d  =  1.0  cm)  shows  that  the  peak  frequency  goes  from 
5.00  hz  down  to  4.94  hz  as  the  amplitude  goes  from  0.2  cm  up  to 
0.8  cm.  Fig.  21  (d  11  1.5  cm)  shows  that  the  peak  frequency  goes 
from  5.38  hz  down  to  5.16  hz  as  the  amplitude  goes  from  0.1  cm  up 
to  0.7  cm.  Fig.  22  (d  *  1.7  cm)  shows  that  the  peak  frequency 
goes  from  5.4  hz  down  to  5.20  hz  as  the  amplitude  goes  from  0.2 
cm  up  to  1.0  cm.  Fig.  23  (d  “  2.0  cm)  shows  that  the  peak 
frequency  goes  from  5.43  hz  down  to  5.03  hz  as  the  amplitude  goes 
from  0.2  cm  up  to  1.0  cm  in  steps  of  0.2  cm.  Fig.  24  shows  that 
the  peak  frequency  goes  from  5.61  hz  down  to  5.18  hz  as  the 
amplitude  goes  from  0.2  cm  up  to  0.9  cm. 

Fig  25  is  a  plot  for  the  case  of  d  =  0.8  cm.  Within  the 
experimental  error  the  peak  frequency  remains  constant  at  4.78 
hz,  as  the  amplitude  of  the  wave  increases  from  0.1  cm  to  0.8  cm. 
It  was  experimentally  found  that  at  d  -  0.8  cm  changes  in  depth 
of  0.1  cm  produced  hardly  detectable  shifts  of  peak  frequency. 
It  is  concluded  that  it  is  at  this  depth  (d  *  0.8  ±  0.1  cm)  that 
the  crossover  from  right  leaning  to  left  leaning  occurs.  The 


calculated  value  of  v#  is  4.81  hz  from  Eq.  (  4  ). 

Figure  26  is  a  reduced  plot  of  the  combination  of  Figs.  19, 
23,  and  25.  The  dash  dot  curves  are  theoretical  peak  values  from 
Eq.  (E-24) . 

Fig.  27  contains  plots  of  the  nonlinear  coefficient.  A, 
versus  fluid  depth,  d.  The  top  curve  is  for  the  case  without 
surface  tension  correction  i.e.  Eq.  (15)  and  the  bottom  curve 
is  for  the  case  with  a  surface  tension  correction  i.e.  Eq. 
(16).  As  shown  by  Eq.  (14)  the  nonlinear  coefficient.  A,  for 
various  depths  of  liquid  can  be  calculated  when  the  measured 
values  of  v0,  v0(£)  and  £  are  plugged  into  the  Eq.  (E-22) .  Thus 
the  experimental  value  of  A  can  obtained  from  the  data  of  Figs. 
19  -  25,  and  represented  by  the  points  on  Fig.  27.  From  the 
curves  we  can  get  an  independent  determination  of  the  crossover 
depth  to  compare  with  the  above-mentioned  experimental  value  (0.8 
cm).  This  is  the  value  of  d  at  which  A  vanishes.  The 
intersection  of  a  horizontal  line,  A  *  0,  with  the  top  curve 
yields  a  value  of  d  *  0.85  cm  and  for  the  bottom  curve,  d  =  1.03 
cm. 

From  Figs.  26  and  27,  we  can  see  that  the  curve  for  Eq.(16) 
is  closer  to  the  experimental  data  than  the  curve  for  Eq.(15)  at 
deep-water  region.  The  reverse  is  true  for  the  crossover  depth. 
Overall  both  equations  work  better  in  the  deep-water  region  than 
shallow-water  region.  The  reason  for  this  is,  as  pointed  by  A. 


Larraza,  a  result  of  a  breakdown  of  the  pertubatlon  theory  as  the 
system  approaches  the  nondispersive  liait  at  zero  depth. 

In  summary,  experimental  measurements  shows  that  a 
rectangular  section  resonator  (V  =  2.6  cm)  has  an  amplitude 

dependent  tuning  curves  with  a  peak  frequency,  v  ,  such  that 

/dC  <  0,  for  d  >  0.8  cm, 

P 

dv„/dC  =  0,  for  d  =  0.8  ca, 

P 

dvVdC  >  0,  for  d  <  0.8  ca. 

P 

The  nonlinear  coefficient.  A,  as  a  function  of  depth,  d,  is  found 

to  be  in  reasonable  agreement  with  the  theoretical  results  of 

(2)  (j) 

Miles  ,  and  Larraza  and  Putterman'  .  As  to  the  discrepancy 

between  theory  and  experiment,  especially  at  the  high  amplitude 
liait.  Professor  Putteraan  expressed  his  satisfaction  with  the 
extent  of  agreement.  He  often  emphasizes  the  theory  is  a 
perturbation  calculation  good  only  to  third  order  in  <<  1, 
which  is  much  smaller  than  the  amplitudes  used  in  the  experiment. 


CHAPTER  V 


THE  FREE  DECAY  RATE  OF  THE  NONPROPAGATING  SOLITON 

A.  WATER- SURFACE- WAVE  DAMPING  IN  CLOSED  BASINS 

The  damping  of  water-surf ace-wave  in  a  closed  basin  has  been 

extensively  studied  by  Benjamin*1***,  Ursell*12*,  Case, 

Parkinson*18*,  Keulegan*19* ,  Van  Dorn*20*  and  Miles*21* .  If  we 

nse  6^  and  &B  to  respectively  represent  the  decay  rate  due  to  the 

viscons  disippation  in  the  neighborhood  of  the  rigid  wall  and  in 

the  neighborhood  of  the  free  surface,  and  &L  to  denote  the  decay 

rate  due  to  the  capillary  hysteresis  associated  with  the  meniscus 

formed  at  the  contact  line  between  the  free  surfaee  and  the  rigid 

wall,  in  general  6^  will  be  much  smaller  than  the  other  two  for 

our  case.  It  was  shown  theoretically  by  Ursell*1***,  Case  and 
(18) 

Parkinson'  that  for  a  clean  bounded  surface  &a  is  negligible 
comparing  with  While  experimental  measured  value  of  may 

be  as  large  as  three.  Van  Dorn'  '  clarified  the  discrepancy  by 
his  experiment,  he  found  that  the  decay  rate  agreed  with  6^  when 
the  water  was  fresh,  then  the  former  tended  to  increase  with  time 
to  some  higher  limiting  value  within  an  hour  or  so.  He 
attributed  the  observed  higher  value  of  to  a  surface  film 
produced  by  spontaneous  contamination.  Considering  these  effects 


Miles'^'1'  derived  the  following  eqnetions: 


■  2J“w)7(1,5 


8,  -  jk^)*  (C  -C.  )coth(^) 
e  4  pu  r  i  V 


(  17  ) 


(  18  ) 


8l  «  ^^(l+^JIn/tpgLi;)] 


(  19  ) 


where  ji  is  the  surface  tension,  L  is  the  length  of  the  basin,  q 
is  the  viscosity,  (C^-C^)  the  ratio  of  the  surface-film  damping 
to  that  which  would  he  produced  by  an  ineztensible  film,  may  be 
approximately  taken  between  0  and  1  for  low  frequency 
water— surf ace-wave  depending  the  extent  of  contamination  at  the 
surface,  k  is  a  semiemperieal  coefficient  which  is  related  to 
capillary  hysteresis  at  the  meniscus  on  the  wall,  it  may  be 
approximated  by  0.05  for  a  hydrophilic  basin  (e.g.  clean  water 
on  glass) . 
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A  sealed  glass  channel,  filled  with  100%  ethyl  alcohol  np  to 
2  cm,  with  the  inside  dimensions  of  L  38x  V  2.6x  E  7.0  in 
centimeters  is  nsed  to  measure  the  free  decay  rates  of  the  low 
amplitude  oscillation  of  the  (0,1)  mode  as  well  as  the  solitons. 

Let  us  first  to  calculate  the  decay  rate  of  the  (0,1)  mode 

by  using  the  abowe-mentioned  formula  worked  out  by  Miles. 

Substituting  the  following  data  into  the  Eqs.  (  17  )  to  (  19  ): 

L  *  38cm,  V  ■  2.6cm,  d  “  2.0cm 

* 

p  •»  23  dyne/cm,  p  =  0.79  g/em  ,  q  -  0.012  g/(sec*cm) 
k  -  0.05,  u  «  2nx5.1  rad/sec,  0<(Cj-Ci)<l 

we  get 

6f  *  0.2  sec  1 
0  <  6g  <  0.3  sec-1. 

6^  "  0.02  sec  1. 

Thus  0.2  sec-1<  8  <  0.5  sec-1. 

Since  we  known  from  experiments  that  the  main  dissipation 
comes  from  the  walls  and  the  surface  of  the  liquid,  the  channel 
is  carefully  cleaned  by  using  the  same  procedure  as  the  one  we 
used  to  clean  the  small  basin  before  the  tuning  curves  were 


measured.  Sight  after  being  filled  vith  100%  ethyle  alcohol  as 
the  working  fluid,  the  channel  it  sealed,  thus  the  contamination 
at  the  liquid  surface  would  be  minimized. 

Fig.  28  is  the  plot  of  a  decaying  (0,1)  mode  surface-wave 
with  the  frequency  of  5.1  hz  and  the  initial  amplitude  of  0.5  cm. 
Decay  rate  is  determined  to  be  0.4  sec”*,  which  is  in  the  range 
of  the  theoretical  prediction. 

«e 

The  free  decay  rate  of  the  soliton  was  determined  in  two 
different  ways.  In  one  the  channel  is  driven  by  the  oscillating 
table  in  its  width  direction,  and  once  the  soliton  is  created  and 
reaches  a  stable  state,  the  driving  source  is  turned  off.  The 
free  decaying  oscillation  data  were  stored  in  a  wave  analyser 
(6000,  Data  Precision)  and  plotted  out  by  a  z-y  recorder.  In  the 
other  method  the  channel  is  oscillated  vertically  at  twice  the 
soliton  frequency  by  a  vibration  exciter  (4809,  BrueldKjaer) . 
Identical  decay  rates  were  obtained  for  both  methods. 

Fig.  29  is  a  plot  of  the  free  decay  soliton  with  the 
initial  peak  amplitude  of  1.8  cm  and  the  frequency  of  5.1  hz. 
Comparing  with  Fig.  28  clearly  the  soliton  decays  faster  than 
the  corresponding  (0,1)  mode  initially.  Fig.  30  is  a  plot  of 
the  ampltude  versus  time.  The  points  are  the  amplitudes  of  the 
same  soliton  as  Fig.  29,  and  the  circles  are  the  amplitudes  of 
the  same  (0,1)  mode  surface-wave  as  Fig.  28.  From  the  plot  we 
can  see  the  initial  decay  rate  of  the  soliton  (frequency  “  5.1 


hx,  aaplitude  ■  1.8  ca)  is  0.7  sec  *  and  changes  to  0.4  sec”*  as 
its  amplitude  drops  below  0.S  cm,  which  is  seen  to  be  the  saae  as 
the  decay  rate  of  the  corresponding  saall  aaplitnde  transverse 
■ode.  The  possible  expanation  of  the  higher  initial  decay  rate 
is  that  the  soliton  has  very  rich  higher  haraonic  coaponents 
which  in  general  decay  faster  than  the  fundaaental  of  the  (0.1) 
node. 
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CONCLUSION 


'  'The  lavs  of  physios 

should  be  staple. '' 
A.  Einstein<22) 


A.  SUMMARY 

The  experimental  discovery  of  a  unique  nonpropagating 
hydrodynasiic  soliton  whose  excitation  is  that  of  a  transverse 
surface-wave  in  a  channel  is  reported.  Like  the  polaron  in 
solids,  its  energy  is  highly  localized  through  a  self-trapping 
effect.  The  profile  of  the  soliton  can  be  decribed  by  a 
hyperbolic  secant  function  and  the  self-trapping  occurs  because 
the  cutoff  frequency  for  the  transverse  node  is  a  decreasing 
function  of  the  wave  amplitude.  Accordingly  the  appropriate 
frequency  for  generation  of  the  peak  energy  at  the  center  of  the 
soliton  is  below  the  cutoff  frequency  at  the  evanescent  wings  and 
above  the  aaplitude  dependent  cutoff  frequency  at  the  peak. 
Measureaents  also  show  that  this  self-trapping  behavior  is  absent 
when  the  liquid  depth  falls  below  a  certain  value,  which  is  in 
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disturbance  cannot  escape  bnt  exponentially  decays  into  the 
surrounding  region  foraing  the  soliton.  One  possible  candidate 
Bight  be  the  flexure  node  in  thin  plates. 

Another  possibility  is  to  look  for  the  nonpropagating 
soliton  in  an  acoustic.  optical,  or  even  biological 
waveguide-like  systea.  Of  couse  the  first  thing  to  check  is  to 
see  if  the  systea  has  the  properties  which  satisfy  the 
above-aentioned  conditions. 

There  are  as  yet  no  theoretical  explanations  of  the 
breather-like  oscillations  of  the  in-phase  soliton  pair  and  the 
initial  free  decay  rate  of  the  soliton. 


APPENDIX  A 


J.  SCOTT  RUSSEL'S  REPOST  ON  HIS  DISCOVERT 

|  J.  Scott  Russel's  Report  on  Waves  st  the  British 

Association  for  the  Advancement  of  Science  was  first  published  in 
1845^.  The  following  paragraph  is  the  most  fascinating  part  of 
j  the  report  of  his  great  disvovery: 

''I  believe  I  shall  best  introduce  this  phenomenon  by 
describing  the  circumstances  of  my  own  first  acquaintance  with 

r 

it.  I  was  observing  the  motion  of  a  boat  which  was  rapidly  drawn 
along  a  narrow  channel  by  a  pair  of  horses,  when  the  boat 
suddenly  stopped  -  not  so  the  mass  of  water  in  the  channel  which 
it  had  put  in  motion;  it  accumulated  round  the  prow  of  the  vessel 
in  a  state  of  violent  agitation,  then  suddenly  leaving  it  behind, 
rolled  forward  with  great  vilocity,  assuming  the  form  of  a  large 

» 

solitary  elevation,  a  rounded,  smooth  and  well-defined  heap  of 
water,  which  continued  its  course  along  the  channel  apparently 
without  change  of  form  or  diminution  of  speed.  I  followed  it  on 
horseback,  and  overtook  it  still  rolling  on  at  a  rate  of  some 
eight  or  nine  miles  an  hour,  preserving  its  original  figure  some 
thirty  feet  long  and  a  foot  and  a  half  in  height.  Its  height 
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gradually  diminished,  and  after  a  chase  of  one  or  two  miles  I 
lost  it  in  the  windings  of  the  channel.  Such,  in  the  month  of 
August  1834,  was  my  first  chance  interview  with  that  singular  and 
beautiful  phenomenon  which  I  have  called  the  Wave  of 


APPENDIX  B 


KdV  TYPE  SOLITONS 


In  1895  Korteweg  and  da  Vrias  inclndad  both  dispersion  and 
nonlinearity,  and  derived  the  ao-oalled  KdV  equation.  There  are 
different  versions  of  derivation  for  the  equation.  The  one, 
introduced  as  follows,  is  worked  out  by  K.  Olsen,  H.  smith  and 
A.  C.  Scott  in  1984 <23). 

For  an  incompressible,  irrotational  and  invisoid  fluid  in  a 
water  tank,  the  velocity  potential  (Fig.  3),  0,  satisfies  the 
Laplace  equation: 


V*  0  «  0 


(B-l) 


Let  us  consider  waves  propagating  in  z  direction  in  the 
water  v.«nk  of  depth  of  d.  The  displacement  of  the  free  surface 
is  specified  by  the  function  £(x,t)  defined  by 


Boundary  condition  at  the  bottoa  of  the  tank  it  the  velocity,  v, 
must  be  zero  in  z  direction: 

(v  )  .  =  (&0/6z)  ,  =  0  (B-3) 

x  x  *  -d  x  *  — d 

At  the  free  surface  the  boundary  equations  are  specified  the 
following  two  equatioins: 

&0/&z  —  80/6x*45/x  ■  =  0  (B-4) 

»0/bt  +  j(v0)*  +  gC  =  0  (B-5) 

where  the  first  one  is  the  kinematic  boundary  equation  and  the 
second  is  Bennoulli's  equation. 

Let  us  express  a  solution  of  Eq.  < B— 1 )  satisfying  the 
boundary  conditions  of  Eq.  (B-3)  as  a  wave  packet: 

0  »  2^/^dkcoshk(d+z)exp(ikx)g(k, t)  (B-6) 

where  the  weight  function  g(k, t)  has  its  fourier  transform  f(x,t) 


equal  to  the  value  of  0  at  the  bottom  of  the  tank,  z  =  -d.  we 
perform  a  simultaneous  expansion  for  z  **  ^(x.t)  of  Eq.  (B-6)  in 
terms  of  o  -  ?MX/d  and  0  *  (d/A)*  -  0(a),  where  A  is  the 

characteristic  length  of  the  soliton,  and  get  that 


0  2  -  ^(d  +  0  fu 


(B-8) 


r  -<<  *  0'u  ♦  J«*  ♦ 


(B-9) 


Thus  Eq.  (B-5)  becomes 


V  f*xx,  -  0 


(B-10) 


by  keeping  terms  to  a  *  ^max^  *n^  ^  onr 

expansion,  while  Bernoulli's  equation  is 


ft  *  «<  -  2d>'ixt  +  K>'  -  ° 


(R-ll) 


From  Eq.  (B-10)  and  Eq.  (B-ll)  the  resulting  linear  dispersion 
relation  for  propagating  plane  wave  becomes 


»  a  a ,  1  *  st 

w  *  c,  It  (1  -  :k  d  ) 


(B-12) 


where  c.  «  zd. 


If  the  function  2  =  f  and  the  dimensionless  valuables  x  = 

7,  t*  =  °  f  and  “*  “  (»c#)/(«CI1|lx)  are  introduced,  Eq. 

'max 

(B-ll)  and  Eq.  (B-12)  can  be  rewritten  as  follows  (dropping  the 
stars) : 


.  IHmOs],  -  -  o 

*t  +  ?x  -  X,t  ♦  «“x  - 


(B-13) 


(B-14) 


In  order  to  determine  an  equation  for  £,  valid  to  first  order  in 
a  and  0,  we  write  2  -  {  +  aA  +  0B  and  insert  in  Eqs.(B-13)  and 
(B-14).  The  two  equations  become  consistent  if  A  *  -  and  B  » 
if 

3  xx,  and  yield  the  KdV  equation: 

«,  * «, +  !•«,  *  K.«  -  (B-15> 

With  the  time  and  space  variables  restored  this  equation  becomes 
the  KdV  equation  in  its  usual  form: 


*t  +  c*?x  +  f(c*/dK4x  +  6d*C#;xxx  "  °*  (B"16) 

It  is  convenient  to  nondimensionalizo  Eq.  (B-16)  in  to  a 
reference  frame  moving  with  the  basic  wave  speed  ov  by 
introducing  the  following  variables: 
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?  *  (x-c0t)/A,  t  =  j(j)  (c,t/A) ,  i\  * 


(B-17) 


The  end  result  is 


HT  +  +  “  °*  (B-18) 

The  soliton  solution  of  Eq.  (B-18)  is 

1 

q  -  C-4xMch*I(;iijix/2)*(5  _2Cmax,c)l  (B-19) 

lor t ewe g  and  de  Vries  derived  a  somewhat  more  general 

equivalent  of  Eq.  (B-16),  in  which  they  included  a  surface 

tension  g,  suoh  that  coefficient  d  in  Eq.  (B-16)  is  replaced  by 

L[l-(^)*/d].  If  d,  <  (^)*  (0.5  cm  for  water),  A  must  he 
PS  i  PS 

replaced  by  |Atl-(^)  /d]  I ,  and  Eq.  (B-18)  becomes: 

\  ~  +  =  °*  (B-20) 

Now  the  solution  for  this  equation  is  a  solitary  wave  of 
depression  rather  than  elevation. 


Fron  Eq.  (B-19)  we  notice  that  the  velocity  of  the  soliton 
is  related  to  its  amplitude  in  such  a  way  that  the  soliton  would 


APPENDIX  C 


F.P.D.  RECURRENCE 


After  world  war  II,  Feral  be cane  interested  in  the 
developnent  end  potentialities  of  the  electronic  coapnter.  The 
first  physical  problem  he  decided  to  try  is  a  problea  for 
heuristic  work.  An  one  dimensional  dynamic  system  of  64 
particles  with  fixed  boundary  condition  and  nonlinear  coupling 
between  nesst  neighbors  was  studied  on  the  Los  Alamos  computer 
Maniac  I.  The  nonlinear  terms  considered  are  quadratic,  cubic 
and  broken  linear  types.  The  results  are  analyzed  into  fourier 
components  and  plotted  out  as  a  function  of  time.  The  original 
motivation  was  to  observe  the  rates  of  thermal ization.  The 
results  of  the  calculations  were  quite  surprising;  the  time 
evolution  of  energy  shows  very  little  tendency  toward 
equipartition  of  energy  among  the  degrees  of  freedom. 

Fig.  32  is  a  plot  of  the  energy  of  the  system  versus  time. 
The  nonlinear  coupling  force  here  is  quadratic.  The  initial 
position  of  the  64  particles  string  is  a  pure  sine  wave  (the 
first  mode).  At  beginning,  the  gradual  energy  flow  from  the 
first  to  the  second  and  the  third  etc.  was  observed.  Later  on. 
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however,  this  gradual  sharing  of  energy  among  successive  nodes 
ceases.  Instead,  it  is  one  or  the  other  node  that  predoninates. 
It  is  only  the  first  few  nodes  which  exchange  energy  among 
themselves.  Finally  mode  I  cones  back  to  within  few  persent  of 
the  initial  value.  This  kind  recurrence  occurs  to  the  cubic  and 
broken  linear  coupling  too.  In  1970  M.  Toda  proved  that  this 

kind  phenomenon  can  be  approximately  represented  by  the  soliton 

(24) 

solution  of  the  so-called  Toda  lattice 


APPENDIX  D 


ENVELOPE  SOLITONS 


It  was  theoretically  discovered  by  Lighthill'  in  1965 

that  weakly  nonlinear  (kC  _  »  «  «  1,  k  ii  wave  nunber,  _  is 

the  characteristic  anplitude  of  the  disturbance)  deep-water  (kd 
>>  1,  d  is  the  depth  of  water)  wavetrain  are  unstable  to 

Modulation  pertubations.  In  1968  Zakharov^^  found  that  this 
Modulation  instability  leads  to  the  fact  that  the  tine  evolution 
of  weakly  nonlinear  deep-water  wave  obeys  the  nonlinear 
Schrodinger  equation: 


ij£  +  Ad*A/3q*  *  v |A* lA  (D-l) 

1 

where  n  *  ekCj-c^t),  x  =  s*(gk)*t1  cg  *  f*.  x 
A  «  k  /[2(gk)  ]*dCj/d«,  v  -  ^(^)  -2gh/[k(kh)  (gh  -  cg  )] 

and  the  free  surface  of  water  is 


z  ■  Re([ei(kh)  A/k  ]ezp[i(kx-<i>t].  (D-2) 

(26k 

In  1971  Zakharov  and  Shabat  solved  NLS  exactly  by  using 

(27) 

the  inverse  scattering  technique  developed  by  Gardner  et  al 


The  soliton  eolation  of  Eq.  (D-2)  is 
x 

A  *  al~|  sech[a(q  -2bT]exp[ibT\  +  iA(e  -b  )v]  (D-3) 

The  eolation  shows  that  the  disturbance,  which  is  mowing  with 
group  velocity  c^,  is  a  localized  wave  packet  with  a  fast 
oscillation  of  carrier  frequency  w  and  an  envelope  specified  by 
hyperbolic  secant  function.  So  it  is  known  as  the  envelope 
soliton. 

This  theoretical  prediction  was  verified  by  Tuen  and  Lake's 
experiments in  1975.  They  use  a  0.915  m  x  0.914  m  x  12.9  m 
water  tank  with  a  wavemaker  at  one  end  and  a  wave-absorbing  beach 
at  the  other.  The  wavemaker  is  a  hinged  paddle  activated  by  a 
hydraulic  cylider.  The  time  evolution  of  each  initial  wave  was 
measured  by  capacitance  wave  amplitude  gauges,  which  were  put  at 
1.53,  3.05,  4.58,  6.10,  7.63  and  9.15  m  in  downstream  of  the 

wavemaker.  Experimental  results  show  that  any  initial  wave 
packet,  other  than  the  soliton  envelope  described  by  Eq.  (D-3), 
eventually  evolves  into  a  number  of  envelope  solitons 
Measurements  also  show  solitons  survive  from  the  interaction  of 
two  envelope  solitons  (Figs.  33  and  34). 


FIGURE  34  The  interaction  of  envelope  solitons.  Left-hand 
trace:  first  pnlse  alone.  Center  trace:  second  pulse  alone. 
Right-hand  trace:  the  second  one  overtakes  and  passes  through 


the  first  one. 
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APPENDIX  E 


THE  THEORY  OF  NONPROPAGATING  SOLITONS 

The  theory  of  nonpropa gating  solitons  were  worked  oat  by 

(2)  (3) 

Miles  ,  aad  by  Larraza  aad  Putteraan  independently.  Larraza 

and  Pnttenan's  approach  neglects  friction  and  drive,  while 

Miles's  approach  inclndes  dissipation  as  well  as  paraaetric 

driwe,  essentially  identical  resnlta  were  achieved.  This  fact 

iaplies  that  neither  dissipation  nor  paraaetric  driwe  are  central 

to  onr  problea.  For  the  sake  of  siaplieity  only  the  firat 

approach  is  introduced  as  follows. 

For  an  inooapressible  inwiscid  and  irrotational  flnid  in  a 
gravitational  field,  if  surface  tension  effects  are  neglected  the 
velocity  potential,  #,  satisfies  Laplaee's  equation. 

▼  0  -  0.  (B-l) 

For  the  problea  we  want  to  solve  (Fig.  3)  the  following 
boundary  condition  should  be  satisfied: 


At  the  free  surface  x  *  C(x.y.t)  the  boundary  conditions  are 
specified  by  two  equations.  One  is  so-oalled  kinenatic  boundary 
condition,  which  is  due  to  the  fact  that  wave  notion  be cones 
stationary  when  it  is  observed  in  a  reference  frane  noving  with 
it.  Consequently 


d/dt[z-{(x,y,t)]  ■  0  (E-4) 

or 

Ct+*xCx+0y<;y  -  0X.  »t  x  -  C(x.y.t) . 

The  other  is  the  d'naic  boundary  condition  which  is  specified  by 
Bernoulli's  equation: 

*t+»5+l/2(*0)*  -  0,  at  s  -  {(x.y.t).  (E-5) 

If  the  surface  displacenent  {(x.y.t)  is  elininated  in  favor 


of  the  velocity  potential,  0,  the  following  equation  for  0  valid 
up  to  terns  that  are  cubic  in  derivatives  of  0  can  be  derived: 


"it***.  "  -{<*»>  Y,<v#t«/*)v,l 

-1/2{<W0)  *^+(00)  *y*y+<*#>  >  S% 

-2/*[(V0)‘t0t]x  -{0x+0tt/g)1[(70)‘ 

-2/g(0t*)xl)  +0(*4)»  at  x  ■  0. 

Up  to  tens  that  are  quadratic  ia  gradients  of  0,  Eq.  (E-6)  can 
be  rewritten  as 


(E-6) 


g(  ■  “0t+l/2[l/g(0t)1x-(w0)*l  at  *  -  0.  (E-7) 

Considering  the  weakly  nonlinear  problen  of  a  disturbance  with  a 
high  frequency  of  notion  «  in  y-direetion  nodulated  by  an 
enr elope  Cx(x, t)  ia  s-directioa  and  satisfying  the  following 
nultiple-seale  requirements: * 

-  •  «  1,  CB-«) 

(dlog(/dx)/(dlog{/dy)  ■  0(s),  (E-9) 

-  0(«*) ,  (E-10) 

[i/(w;1).a;1]/3t  -  o(s*),  (e-ii> 

where 
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* nJ\  ».w  *J  xAa 


\ 

V*.-* 


Mx  -  gkT,  T  -  tanh(kd),  k  -  n/W.  (E-12) 

The  solution  of  Laplace's  equation  specified  by  tke  boundary 
conditions  Eqs.  (E-2)  and  (E-3)  is  given  by 

0s  (0x(x, t) [coshk(z+d)/coshkd]«cosky*exp(i<*t)+c.c. )  (E-13) 

+  l0#  ^  (x,  t)exp(2it»t)+c.c.  ]+  [cosh2k(z+d)/eosh2kd]«cos2ky 
•(I01(x,t)exp(2i**t)+c.c.]+0s<*)(x.t)j  ♦<C-l/(2k)]ed*01/ax* 

•( cosky /cosbkd)*exp(iwt) [xsinhk(z+d)-dexp(-k(z+d) )]+c.c.} 

+0# (x+i ( x+d) , t J+0, [x-i( z+d) « t}+0( e4 ) . 

Substituting  Eq.  (E-13)  into  Bq.  (E-7)  and  equating  to  zero 
coefficients  of  the  sane  cos(nky)exp(niwt)  dependence,  the 
following  solutiona  of  different  orders  can  be  detemined 

■  kgT  (valid  to  0(e))  (E-14) 

0,  -  [(-3ik,01*)/(8-Ta)](l-T4)  (valid  to  0(s*)>  (E-15) 

0,(>)  -  Iik*/(8«)  10/(1+31*)  (E-16) 

(valid  to  0(e*)) 

0.(#)  -  0  (valid  to  0(s*)) 


(E-17) 


(E-18) 


2i»d0l/dt-c  9  01/dz  +(*i 
-*»%)0i-Ak  **  *x  -  0  ( valid  to  0(«*)) 

where 

c*  -  jjlT+kd ( 1-T* ) ]  (E-19) 

end 

A  -  l/8(tfT4-5T*+16-9T"*).  (E-20) 

If  the  surface  tension  effeets  are  included,  as  worked  ont 
by  Miles Eqs.(E-19)  and  (E-20)  should  respectively  Modified 

as 


•.*  -  (kg  +  (E-21) 

P 

and 

A  -  |[  1+ ( 1-T* ) *+J( l+4(i* )_1(1+T*)*  -  J(T*-  4(i*)_1(3-T*)*J  (E-22) 

where  p  is  surface  tension,  and  p*  *  p(j|)  /(pg). 

Ve  introduce  aaplitude  dependent  cutoff  frequency,  wt(£). 


which  is  given  by 


V.<W**.(1-AI  /<128f  v«  >•  <E“24> 

Under  the  condition  that  A>0,  which  neens  that  d>1.022/k  froa 
Eq.  (E-20) ,  •,<;)<»*..  Aa  ».<w<«,K),  Eq.(E-18)  possesses  a 

aoliton  aolntion  of  the  fora 

i 

*x  "  1 2 ( »# *-*•* ) / (Ak4 ) ] *sech [ (  a/Vlx/o1]  (E-25) 

The  anrface  diaplacenent,  ^(x,y,t),  can  be  calcnlated  froa 
Eq. (E— 7)  aa 

((x,y,t)  -  ^([-ia#xexp(ittt)eosky+e.e.]  (E-26) 

+[Jk*l*x|*<T%l)cos2ky  ♦  Jk*l0xl*(T*-l) 

”  ^k  #x  ( j*-l)exp(2i«t)coa2ky+c. c. ] } 


In  suaaary,the  aolntion  haa  the  following  featnrea: 

(1)  the  ahape  of  the  modulation  is  giwen  by  hyperbolic  secant 
function. 


(2)  the  aoliton  is  stationary  along  x-direetion 
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